function ds_aug = augment_dataset(ds)
ds_aug = ds;
for i=1:size(ds,2)
   for j=i:size(ds,2)
      ds_aug = [ds_aug ds(:,i).*ds(:,j)]; 
   end
end
end